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Abstract — In this paper, we discuss iterative channel estimation 
for the Long Term Evolution (LTE) downlink. The LTE standard 
utilizes coherent detection and therefore to achieve a high data 
rate over mobile radio channels, it is essential to obtain high 
accurate channel state information at the receiver side. For 
channel estimation purposes, the LTE standard provides training 
data known as pilot symbols. The open question is whether 
the accuracy of channel estimation based on pilot symbols is 
sufficient to achieve such high data rate transmission. The quality 
of the channel estimates can be further enhanced if soft or hard 
information bits, which are available after the decoding stage, 
are utilized by the channel estimator to re-estimate the channel. 
In order to quantify the gain provided by the iterative channel 
estimation, several scenarios have been considered in this work. 

Index Terms — LTE, OFDM, Channel Estimation, Iterative 
Receiver. 

I. Introduction 

Channel estimation plays a crucial role in the performance 
of wireless communication systems, since its knowledge is 
utilized to detect the data symbols [1]. With enhancing the 
channel estimation performance, also the performance of the 
entire system is improved [2]. Many investigations have been 
conducted on channel estimation in LTE [3], [4], [5]. The 
majority of them investigate non-iterative channel estimation 
algorithms based only on the known pilot symbols provided by 
the LTE standard. In recent years, receivers with iterative chan- 
nel estimation have become more and more attractive because 
of their superb performance [6]. Hard or soft information bits 
from the decoding stage are utilized to estimate the data sym- 
bols and then applied to further improve the channel estimation 
[7], [8]. In this paper we discuss and analyse the performance 
of different iterative channel estimation algorithms for the 
block fading case. 

The main contributions of the paper are: 

• We calculate the soft estimated data symbols for the LTE 
standard (according to [9]) for 4-Quadrature Amplitude 
Modulation (QAM), 16-QAM and 64-QAM constella- 
tions symbols. 

• We derive various iterative channel estimation algorithms 
with different complexities and performance behaviours. 

• We investigate the impact of processing the soft (a- 
posteriori or extrinsic) and hard feedback information for 
channel estimation. 

The remainder of the paper is organized as follows. In 
Section II, the system model is described. Various iterative 
channel estimation algorithms are presented in Section III. 



In Section IV, we evaluate the performance of the channel 
estimators in terms of data throughput and Mean Square Error 
(MSE). Finally, Section V concludes the paper with a summary 
of most important points obtained in this work. 

II. System Model 

In this section, we briefly give an overview of the pilot 
symbol structure defined in the LTE standard and discuss the 
structure of an iterative receiver. Afterwards, we introduce a 
signal model and we show how the soft symbols are calculated. 

A. Overview 

In LTE systems, the initial channel estimates can be ob- 
tained based on pilot symbols, which are placed in the 
Orthogonal Frequency Division Multiplexing (OFDM) time- 
frequency grid. In Figure 1, the pilot symbol structure for one 
subframe corresponding to normal cyclic prefix configuration 
for up to four transmit antenna ports is shown (for differ- 
ent antenna ports different color boxes) [9]. Pilot symbols 
transmitted from multiple antenna ports are orthogonal to 
each other. That is, when one antenna port transmits a pilot 
symbol, the other antenna ports do not transmit. For the 
iterative channel estimation we take advantage of hard or soft 
estimation of data symbols (white boxes), available after the 
decoding stage, which in addition to the pilot symbols are 
utilized to re-estimate the channel after the initial estimation. 
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Fig. 1. Pilot symbol structure. 

B. Iterative Receiver 

Figure 2 depicts the iterative receiver structure. The channel 
state information knowledge, derived from the pilot symbols, 
is utilized for demodulation and soft de-mapping of the OFDM 
symbols. Afterwards, the detector provides to the decoding 
unit soft information in terms of Log Likelihood Ratios (LLRs) 
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for each transmitted information coded bit. In order to take 
advantage of the feedback information, the decoding unit has 
to be modified and new components such as the soft symbol 
mapper has to be added to the receiver (the iterative receiver 
part in Figure 2). The modified decoding unit, besides hard 
decoded bits, outputs also soft estimated coded bits. This 
feedback information is utilized by the (soft) symbol mapper 
to calculate the estimated data symbols, which in turn, together 
with the pilot symbols are applied by the channel estimator to 
re-estimate the channel. The process is performed again and 
again and new channel estimates are obtained. Through the 
successive iterations, the overall system performance can be 
improved. 
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Fig. 2. Iterative receiver structure. 
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C. Signal Model 

The received frequency domain signal r„ r „ B (after OFDM 
demodulation) at the receive antenna port n Y and OFDM 
symbol n s is obtained as 



E« 

n t = l 



n t ,n e h ratjIlr 



(1) 



where S„ t ,„ B = diag(s 

n t ,n s ) is a diagonal matrix and s„ t „ B 
is the transmit signal vector at the transmit antenna port rit 
(n t = 1 . . . Nt) and OFDM symbol ?i s , which comprises the 
pilots s PtntinB>nk and data subcarriers Sd, nt ,n s ,n k , where n k is 



the subcarrier index. The vector h,- 



is the channel vector 



in the frequency domain and n„ r „ s is the additive zero mean 
white Gaussian noise vector with the variance a^- The length 
of vectors s„ tj „ B , h„ t „ r and n„ r .„ B is Nk, which corresponds 
to the number of subcarriers within an OFDM symbol n s . If we 
concatenate Equation (1) for Ns OFDM symbols, we obtain 



h 



(2) 



where S„ t is of size (Nk Ns) x Nk, and r„ r and n nr are 
of size (Nk Ns) x 1. On the one hand, for initial channel 
estimation the sum in Equation (2) disappears because if 
there is a pilot symbol located at the n t -th transmit antenna 



port within the n s -th OFDM symbol on the n^-th subcarrier, 
symbols at the remaining transmit antenna ports at the same 
positions are zero. Therefore, Equation (2) extracted at pilot 



symbols p is simply r Pt , 



On the 



other hand, for iterative channel estimation, the matrix S„ tj „ B 
from Equation (1) is replaced by the matrix S„ t .„ B , in which 
the data symbols are replaced by the estimated soft or hard 
data symbols. Thus, Equation (2) for the iterative case can be 
written as 



Sh„ 



(3) 



with S, of size (Nk Ns) x (Nk Nt), comprising the pilot 
symbols and the soft or hard estimated data symbols. 

D. Soft Symbol Calculation 

The a-posteriori or extrinsic LLRs of coded bits from the 
output of the decoding unit are considered as a-priori informa- 
tion and utilized by the soft symbol mapper to calculate the 
soft estimated data symbols as [10] 



s d = E [s d ] 



M 

E 

i=l 



Sd; 



N c 



(4) 



where Sdj is one of the constellation points that is formed 
by the information bits /2, • • • , fNa)> M is the number 
of symbols in the symbol alphabet and Nc = log 2 M is the 
number of bits per symbol. P(/„ b (sd.i)) is the probability of 
the bit f Jlb that corresponds to the given symbol constellation 
point Sd.i- For notational simplicity we have omitted the 
transmit antenna port nt, OFDM symbol n s and subcarrier 



iik indexes from Sd ■ 



The LLR of the coded bit f„ b is 



computed by the log-ratio of the probability of the information 
bit being 1 over the probability of the information bit being 0 

P(/»» = l) 



L(/nJ=ln 



We know that P(/„ 



1 



p(/n 6 = oy 

P(/n„ = 0) 



(5) 

1. Using the 



definitions of LLRs above, we find the probability of the 
information bit being 1(0) as 

l(0)XL(/ nf) ) 

P(/n 6 = 1(0)) = J^MKJ- ^ 

Furthermore, Equation (6) can be written in terms of tangent 
hyperbolic functions 

P(/n t =1(0)) = i(l + Rtanh(L(/„J/2)). (7) 

Any Gray-mapped QAM constellation symbol can be divided 
into real and imaginary parts Sd = Re(sd)+jlm(sd) [10]. This 
implies that the real and imaginary components of each soft 
symbol are independent. For the real Re(srf) and imaginary 
lm(sd) part of the soft estimated symbol Sd, by applying 
Equation (7) in Equation (4), for 4-QAM, 16-QAM and 64- 
QAM (for the LTE standard defined modulation alphabet 
according to [9]) we find 
. 4-QAM 



Re(s d ) = -ai and lm(s d ) 



-a 2 . 



(8) 
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. 16-QAM 

Re(Sd) = — aids — 1a\ and lm(Sd) = —0204 — 202. (9) 
. 64-QAM 

Re(Sd) = — 4<zi — 2aia3 — 01030,5 and (10) 

Im(s d ) = -4a 2 - 2a 2 04 - 020406, (11) 

where a nb = tanh (L(/„ 6 )/2) is the tangent hyperbolic func- 
tion as a function of the bit LLR L(f„ b ). The soft estimated 
symbols are then obtained as Sd = Re(.?ci) + jlm(srf). 

III. Iterative Channel Estimation 

In this section, different iterative channel estimation algo- 
rithms considered in this paper are derived. 



A. Iterative Least Square (LS) channel estimator 

Iterative LS channel estimate tries to find a solution h 

(i) 



(0 

ILS,n r 



by minimizing the Euclidean norm squared J(hJ L ; s n ) 



Sh 



(») 

ILS, 



The minimum is found by computing the 



gradient of the cost function V J(h ILS ) and equating it to 
the zero vector. Thus, the iterative LS channel estimator is 
derived as [11] 



(12) 



under the condition that rank(S) = NrNt [12], where i is 
the iteration index number. For notational simplicity, we have 
omitted the iteration index i in S^K 

B. Iterative Linear Minimum Mean Square Error (LMMSE) 
channel estimator 

The iterative LMMSE channel estimator provides estimated 
channel coefficients by minimizing the MSE 



E 



1 ILMMSE,n I 



(13) 



Utilizing the second order statistics of the channel the iterative 
LMMSE estimator is derived, which improves the performance 
of the iterative LS estimator by filtering its estimate as [11] 



ILMMSE,n r 



,(Rh„ 



TLS,n r > 

(14) 

where Rh„ ,h„ is the autocorrelation matrix of the channel 
vector h„ r , which is assumed to be perfectly known. 

C. Iterative Approximate LMMSE (ALMMSE) channel estima- 
tor 

The computational complexity that comes with iterative 
LMMSE is very high because of the matrix inversion in- 
volved in Equation (14). Therefore, its use in a real-time 
implementation is limited and restricted to a low complexity 
implementation. In this subsection, we investigate an approx- 
imate algorithm which reduces the complexity of the iterative 
LMMSE and preserves its performance. This estimator is 
presented in [3] for LTE and [13] for WiMAX but derived 
only for initial channel estimation. Therefore, we will modify 



this estimator for the iterative case. As presented in [3] and 
[13], the main idea of ALMMSE estimator is to approximate 
the correlation matrix Rh.h by utilizing only the correlation of 
the L closest subcarriers instead of using the full correlation 
between all subcarriers as in LMMSE case. Assuming that 
the correlation is frequency independent, the approximated 
L x L correlation matrix R^/h ls obtained by averaging over 
the correlation matrices of sizes L x L [14]. As in [3], first 
the number of L closest subcarriers is defined, over which 
the autocorrelation matrix R^ is approximated. Next, the 



channel vector h 



is defined over the interval of these L 



consecutive subcarriers as 



" , n t ,n r ,l ■ ■ ■ ll n t ,n 
" / n t ,n t ,N K -L + l ' 1 



nt,n r ,k+\— 3—] 
n t ,n t ,Nk 



1 T 



; otherwise 



;k>N K 



L-l 
2 

(15) 

where k is the subcarrier index. The output of the estimator 
for the given interval of L subcarriers becomes = 

u(i).(L)J 



U«,(L),T 
ll N T ,n r 



Furthermore, we derive iterative LS 



hjLs „ over all subcarriers and define it over the interval of 
L consecutive subcarriers as above 

T 



(i),ILS lT 

ll-.L ■ ■ 



H 

fh 



[h 



(t), ILS i T 
N T ,n r \l:L 



( l ),ILS lT 

fe-L ; 

),ILS lT 
l,n, \N K - 



;(i),iLS lT 



L+1:N K 



[h 



(i),ILS lT 

N T ,n 



— j ; fc+r- 

i T 



Nk-L+1:Nk 



r 1 ! 



(16) 



Finally, we calculate the iterative ALMMSE channel estimate 
by filtering the iterative LS estimate h^g „ with the filtering 

matrix A{^ LMMSE as follows hn}'^ 
where the filtering matrix is given by 



a(L) C(0.(l) 
-"-IALMMSE^ILS.n, ' 



,(L) 
IALMMSE 



-r[ l) 



R 



(L) 



(s H s)- 



(L) N 



(17) 

Then, we select the k-th subcarrier for the estimated channel 
vector hk} 



£(i),(L)lT 

k ■ ■ 



[hf WL) ' T 



l n N T ,n r \k 



L+k-N K 



W,(L)lT 
N T ,n v \ \L±±-\ 

^W,(L)lT 



N T ,njL+k-N K 



(18) 



IV. Simulation Results 

In this section, we discuss the performance of the it- 
erative channel estimation algorithms based on simulation 
results. We compare the iterative approaches with the pilot 
based approaches. Naturally, the question whether to apply 
a-posteriori, extrinsic or hard estimated data symbols in the 
feedback loop arises. In the following, we investigate the 
impact of processing this feedback information for channel 
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estimation. We distinguish between initial (pilot based), a- 
posteriori, extrinsic, hard and perfect cases. In the initial case, 
in the plots denoted by It. 1 , the channel is estimated based only 
on pilot symbols. When the system knows the channel state 
information perfectly then this is recognized as the perfect case 
and denoted as PERFECT. The utilization of the a-posteriori, 
extrinsic and hard feedback information in the feedback loop 
is denoted (in the front of the abbreviations of the channel 
estimation algorithms) by app, ext and hard, respectively. 
The channel estimators are compared with each other in terms 
of MSE and throughput versus Signal to Noise Ratio (SNR). 
All results are obtained using the LTE Link Level Simulator 
version "rl089" [15], [16]. Table I presents the most important 
simulation settings. 



Parameter 


Value 


Bandwidth 


1.4 MHz 


Number of subframes 


1000 


Number of iterations 


5 


Number of closest subcarrier L 


12 (for ALMMSE) 


Number of transmit antennas 


2 


Number of receive antennas 


2 


Channel Quality Indicator (CQI) 


4,7,10 


Modulation order 


4-,16-, and 64-QAM 


Transmission mode 


Transmit Diversity (TxD) 


Channel type 


ITU PedB 


Detector 


SSD 



TABLE I 
Simulation settings 

Until otherwise stated, for an iterative ALMMSE estimator 
we chose to use L=12 closest subcarrier for the following 
simulations. Figure 3 depicts the throughput versus the num- 
ber of iterations for the three channel estimation algorithms 
when a-posteriori feedback information is utilized. The highest 
performance gain is achieved from the first to the second 
iteration. Then, from the second iteration on the throughput is 
increased slightly and after the fifth iteration the performance 
is improved negligibly. Therefore, in the following simulations 
only five iterations between the channel estimator and the de- 
coding unit are considered. For further analysis of the iterative 
channel estimation algorithms, we consider only the 4-QAM 
modulation order. The extension to other modulation orders 
is straightforward. For comparison purposes, in Figure 4, the 
throughput versus SNR for 4-QAM, 16-QAM and 64-QAM 
modulation orders, when a-posteriori feedback information is 
utilized, is plotted. As for 4-QAM, approximately the same 
performance pattern are also observed for 16-QAM and 64- 
QAM. 

A. A-posteriori, extrinsic and hard feedback information 

In this subsection, we analyse the performance of the system 
when the a-posteriori, extrinsic and hard feedback information 
are utilized in the feedback loop. Figure 5 depicts the through- 
put (analysed at throughput 0.4 Mbit/s) versus SNR, where 
the channel estimation algorithms (initial and iterative case) 
for the above mentioned feedback information after the fifth 
iteration are compared to each other. We observe that the best 
system performance is achieved when a-posteriori feedback 




3 4 
Number of iterations 

Fig. 3. Throughput versus the number of iterations for a-posteriori feedback 
information. 
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Fig. 4. Throughput for 4-, 16- and 64-QAM modulation orders for a-posteriori 
feedback information. 



information is utilized in the feedback loop. The improvement 
of the a-posteriori case with respect to the extrinsic case is 
approximately 0. 1 dB for the three iterative channel estimation 
algorithms. On the other hand, if we observe the case when 
hard feedback information is utilized we see that the iterative 
hard LS estimator does not improve with respect to the initial 
(pilot based) case, iterative hard ALMMSE estimator obtains 
approximately 0.3 dB SNR gain and iterative hard LMMSE 
estimator 0.15 dB gain compared to their initial cases. Clearly, 
the soft (a-posteriori or extrinsic) case outperforms the hard 
case. A-posteriori LS estimator obtains approximately 0.35 
dB SNR gain, a-posteriori ALMMSE estimator 0.8 dB and 
a-posteriori LMMSE estimator 0.4 dB gain with respect to 
the hard case. Obviously, the major drawback utilizing hard 
feedback information is the error propagation which is more 
severe compared to the soft case. Since the highest system 
performance is achieved when a-posteriori information is 
utilized, in the following, we compare the channel estimation 
algorithms to each other for the a-posteriori feedback infor- 
mation case. 
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B. Comparison of the channel estimation algorithms for a- 
posteriori feedback information 

First, we investigate the ALMMSE estimator in dependence 
of the closest subcarrier index L. In Figure 6 and Figure 7, the 
performance of the iterative (for the fifth iteration) ALMMSE 
estimator for different values of L, and in comparison to LS 
and LMMSE estimators, is plotted. It can be observed that 
with increasing L the MSE decreases while the throughput 
increases. For L = Nk, the iterative ALMMSE estimator 
performance is equal to the iterative LMMSE estimator. It 
is obvious that with increasing L also the complexity of the 
estimator increases. This fact allows adjusting the performance 
and complexity of the estimator to achieve a good trade-off. 
We see that even for L=2 consecutive subcarriers the iterative 
ALMMSE estimator performs better than the iterative LS esti- 
mator. Increasing the number L more than 12 the performance 
of the system is increased negligibly as compared from L=2 
to L=12. Therefore, we conclude that the choice L=12 is a 
good trade-off between complexity and performance. 

Figure 8 depicts the MSE versus SNR for the three iterative 
algorithms (after the fifth iteration) in comparison to their 
initial cases. For the three iterative channel estimation algo- 
rithms, the MSE considerably decreases after five iterations 
compared to the initial case. The MSE decreases most rapidly 
for the LMMSE, followed by the ALMMSE, and then the 
LS. Figure 9 depicts the throughput versus SNR. For iterative 
LS we observe approximately 0.3 dB SNR gain, for iterative 
ALMMSE 1.1 dB and for iterative LMMSE 0.6 dB with 
respect to their initial cases. Iterative ALMMSE outperforms 
the initial LMMSE for about 0.35 dB and loses 0.25 dB 
compared to iterative LMMSE. On the other hand, iterative 
LMMSE channel estimator performs approximately equal to 
the system with perfect channel knowledge, with an SNR loss 
of 0.05 dB. Table II and III summarize in more detail the 
above SNR gain and loss values, measured at throughput of 
0.4 Mbit/s. Each iterative estimator algorithm is compared with 
its initial estimator as well as with the other iterative estimators 
and also with the system with perfect channel knowledge. 




0 5 10 

SNR [dB] 



Fig. 6. Iterative ALMMSE estimator for different values of L compared to 
iterative LS and LMMSE estimators in terms of MSE. 
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Fig. 7. Iterative ALMMSE estimator for different values of L compared to 
iterative LS and LMMSE estimators in terms of throughput. 
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Fig. 8. Comparison of initial (pilot based) and iterative LS, ALMMSE and 
LMMSE estimators in terms of MSE for a-posteriori feedback information. 
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SNR [dB] 

Fig. 9. Comparison of initial (pilot based) and iterative LS, ALMMSE 
and LMMSE estimators in terms of throughput for a-posteriori feedback 
information. 



SNR gain values 


Iter, versus Init. 


Init. LS 


Init. ALMMSE 


Init. LMMSE 


Iter. LS 


0.3 dB 






Iter. ALMMSE 


2.3 dB 


1.1 dB 


0.35 dB 


Iter. LMMSE 


2.5 dB 


1.35 dB 


0.6 dB 



TABLE II 

Summary of SNR gain values of the iterative over the initial 

(PILOT BASED) ESTIMATORS 



SNR loss values 

Iter, versus Iter. | Iter. LS | Iter. ALMMSE | Iter. LMMSE | Perfect 

Iter. TjT - 2 dB ~2.25 dB 2.3 dB 

Iter. ALMMSET - - ~0.25 dB 0.3 dB 

Iter. LMMSE | - - - | 0.05 dB 

TABLE IE 

Summary of SNR loss values of the iterative estimators over 

EACH OTHER AND OVER THE PERFECT CASE 



V. Conclusion 

We showed the calculation procedure of the soft estimated 
data symbols for the LTE standard. Furthermore, we derived 
different iterative channel estimation algorithms with different 
performance behaviours and complexities. We investigated 
the impact of processing the a-posteriori, extrinsic and hard 
feedback information for the channel estimation. The sys- 
tem with iterative channel estimation achieves a considerable 
performance improvement in terms of MSE and throughput 
compared to the system with non-iterative channel estimation. 
The best performance is achieved when a-posteriori feedback 
information is utilized. Among the estimators presented in 
this paper, the iterative LS estimator achieves the lowest per- 
formance. The iterative LMMSE estimator achieves approx- 
imately the performance of the system with perfect channel 
knowledge with an SNR loss of 0.05 dB. Despite the fact that 
the performance of the iterative LMMSE estimator is excellent, 
its complexity is very high for a real-time implementation. 
Therefore, the goal of iterative ALMMSE estimator is to 
reduce the complexity of the iterative LMMSE estimator and 
in turn to preserve its performance. The iterative ALMMSE 
estimator exploits only the correlation between the L closest 
subcarriers. Thus, there is a trade-off between complexity and 
performance by varying the parameter L. Iterative ALMMSE 
estimator with L=12 obtains approximately 1.1 dB SNR gain 
with respect to the initial case and loses approximately 0.25 
dB compared to the iterative LMMSE estimator. 
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